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ABSTRACT. We consider an active scalar equation that is motivated by a model for magneto-geostrophic dy- 
namics and the geodynamo. We prove that the non-diffusive equation is ill-posed in the sense of Hadamard in 
Sobolev spaces. In contrast, the critically diffusive equation is well-posed (cf. 1151 ). In this case we give an 
example of a steady state that is nonlinearly unstable, and hence produces a dynamo effect in the sense of an 
exponentially growing magnetic field. 



The two classical examples of active scalar equations arising fluid dynamics are Burgers' equation from 
compressible fluids, and the two-dimensional vorticity equation from incompressible fluids. The analyt- 
ical study of these seemingly simple equations has generated a substantial amount of new mathematics 
over the past centuries. An active scalar equation from geophysical fluids that has recently received con- 
siderable attention in the mathematical literature (see for instance the reference list in [6]) is the surface 
quasi-geostrophic (SQG) equation. Introduced by Constantin, Majda, and Tabak |[8l as a two-dimensional 
toy model for the three-dimensional fluid equations, the SQG equation has been studied in both the inviscid 
and the viscous versions, and recently it was shown that the critically viscous equation is globally well-posed 
Il3"l l2~l1l . However the possibility of finite time blow-up in the inviscid SQG equation is still open. 

In the present paper we address a class of three-dimensional active scalar equations for which the drift 
velocity is more singular by one derivative than the active scalar. In contrast, the drift velocity for the 
SQG equation is of the same order of derivatives as the active scalar. Our motivation for addressing such 
singular drift equations comes from a model presented by Moffatt l27l for the geodynamo and magneto- 
geostrophic turbulence in the Earth's fluid core. The nonlinear effects in this three-dimensional system 
are incorporated in an evolution equation with singular drift for a scalar buoyancy field. This magneto- 
geostrophic (MG) equation has certain features in common with the SQG equation, most notably that L°° 
is the critical Lebesgue space with respect to the natural scaling for both critically diffusive equations. 
Inspired by the proof of global well-posedness for the critically diffusive SQG equation given by Caffarelli 
and Vasseur 0, we recently used De Giorgi techniques to prove global well-posedness for the critically 
diffusive MG equation [ 15 ]. In contrast, in this present paper we show that the non-diffusive MG and SQG 
equations are distinctly different at the analytical level, not only because of the more singular drift velocity, 
but also because of the structure of the operators relating the velocity to the active scalar. 

We study a class of nonlinear active scalar equations for the unknown scalar field B(x, t), driven by a 
singular velocity vector field U(x,t), namely 



where k > is a physical parameter, and S(x) is a C°°-smooth source term. The velocity U is divergence- 
free, and is obtained from G via 



1. Introduction 



d t O + U ■ V6 

v-jy = o 



S + kAQ 
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for all j G {1, . . . , d}, where {Tij} is a d x d matrix of Calderon-Zygmund operators of convolution type, 
such that didjTijip = for any smooth ip. Here and throughout this paper we use the summation convention 
on repeated indices. For simplicity, we consider the domain to be (x, t) G T d x (0, oo) = [0, 2ir] d x (0, oo). 
Without loss of generality we may assume that Jj d Q(x,t) dx = for all t > 0, since the mean of is 
conserved by the flow. 

The model proposed by Moffatt is derived from the full MHD equations for an incompressible, rotating, 
density stratified, electrically conducting fluid. After a series of approximations relevant to the Earth's fluid 
core, a linear relationship is established between the velocity and magnetic vectors, and the scalar buoyancy 
0. The active scalar equation for 0(x, t) that contains the nonlinear process in Moffatt's model is precisely 
( |1.1| >, but where the divergence-free velocity U is explicitly obtained from the buoyancy as 

Uj = Mj@, (1.4) 

for j G {1, 2, 3}. Here the Mj are Fourier multiplier operators with symbols given explicitly by 

-r> 2^ 2 fe 3 |fc| 2 -(/3 2 /r / )fc 1 fc 2 2 fc 3 n M 

U J 40 2 A||fc| 2 + (/3 2 /r ? ) 2 ^ K J 

M2(fc) = 4^ 2 |fc| 2 + (/3 2 A ? ) 2 ^ ( 6) 

where the Fourier variable k G Z d is such that fe 3 7^ 0. On {A; 3 = 0} we let Mj(k) = 0, since for self- 
consistency of the model we assume that and U have zero vertical mean. It can be directly checked 
that kj • Mf(fc) = for all k e Z d \ {k 3 = 0}, and hence the velocity field U given by ( |1.4| > is indeed 
divergence-free. It is important to note that although the symbols M, are bounded in the region of Fourier 
space where k\ < max{/c2, A; 3 }, this is not the case on the "curved" frequency regions where A; 3 = 0(1) 
and &2 = 0(|/ci| r ), with < r < 1/2. In such regions the symbols are unbounded, since as |fci| — > 00 we 
have 

|Mi(fci, \h\ r , 1)| w l^r, |M 2 (A;i, |fci| r , 1)| « 1^1, |j\? 3 (fci, |fcj| r , 1)| « |fc!| 2r . (1.8) 

In fact, it may be shown that |M(fc)| < C\k\, where C(/3, rj, O) > is a fixed constant, and that this bound 
is sharp. In lTT5ll we make precise the fact that ( |1.1[ )-( [L2"1 ), with velocity given by ( |1.4| )-( fT77| ), is an example 
of the abstract system ( |l.l| )-( fT3| ), by letting be the zero-order fourier multipliers 

Tij = -dii-A^Mj. (1.9) 

We refer to the evolution equation ( |1.1| )-( [L~2] ), with drift velocity U given by ( |1.4| )-( fT77] ), as the magneto- 
geostrophic equation (MG). On the other hand, we refer to the general case when U is given by \\3\ as 
the abstract magneto- geostrophic equation (AMG). We make this distinction to emphasize that some of 
the theorems stated in this paper do not make use of the explicit structure of the symbol of M, whereas 
other results, including ill-posedness, make explicit use of the specific structure given by ( |1.5| >-( [T77| ). In 
addition, we shall refer to the diffusive evolution ( |1.1| ), i.e. k > 0, as MG K (respectively AMG K ), and to the 
non-diffusive case k = as MGo (respectively AMGo). 

The physical parameters of the geodynamo model are the following: O is the rotation rate of the Earth, 
r\ is the magnetic diffusivity of the fluid core, and f} is the strength of the steady, uniform mean part of the 
magnetic field in the fluid core. The perturbation magnetic field vector b(x, t) is computed from 0(x, t) via 
the operator 

bj = (p/ v )(-A)- 1 d 2 M j e, for all j G {1,2,3}. (1.10) 

In the Earth's fluid core the value of the diffusivity n is very small. Hence it is relevant to address both the 
diffusive evolution ( |1.1[ ), and the non-diffusive version where k = 0. As we will demonstrate, the diffusive 
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and the non-diffusive systems have contrasting properties: if k > the equation is globally well -posed and 
the solutions are C°° smooth for positive time |[T5l[T6ll . whereas for k = we shall prove that the equation is 
ill-posed in the sense of Hadamard in Sobolev spaces, but locally well-posed is spaces of analytic functions. 

We first address the non-diffusive version of ( |1.1| >. We prove that the AMGo equations are locally well- 
posed in the space of real-analytic functions. The proof of analytic well-posedness, a Cauchy-Kowalewskaya- 
type result, is given in terms of the Gevrey-class norms introduced by Foias and Temam ifTTTl for the Navier- 
Stokes equations, combined with ideas from Il23ll24l for the Euler equations. Moreover, we show that the 
breakdown of the real-analytic solutions is fully characterized by a Sobolev norm of the solution, and the 
initial data, providing an effective self-consistency tool for numerical simulations of the AMGo equations. 
In addition, we point out that if the matrix Tij is self-adjoint, the operator 6 i->- U is anti-symmetric, and the 
arguments of [5 ] may be used to obtain the local well-posedness of the AMGo equations in Sobolev spaces. 
We emphasize that the specific MGo operator M defined via ( |1.5| )-( fT77] ) is not anti- symmetric. 

In the theory of differential equations, it is classical to call a Cauchy problem well-posed, in the sense 
of Hadamard, if given any initial data in a functional space X, the problem has a unique solution in 
L°°(0, T;X), with T depending only on the X-norm of the initial data, and moreover the solution map 
Y h-> L°°(0, T; X) satisfies strong continuity properties, e.g. it is uniformly continuous, Lipschitz, or even 
C°° smooth, for a sufficiently nice space Y C X. If one of these properties fail, the Cauchy problem 
is called ill-posed. Depending on the specific equation, and on the regularity of the space X considered, 
one or more properties of a well-posed problem may fail, either locally or globally in time. At one end 
of this spectrum we have the dramatic examples of blow-up in finite time for solutions (e.g. iPTOl for the 
Prandtl equations), while at the other end we have the case in which the solution map is continuous with 
respect to perturbations in the initial data, but not C k -smooth for some k > 0, or at least not C°° -smooth 
(cf. |[T7l[T8l[T9ll20l[29i ). In-between cases of ill-posedness may be due to the fact that X is too large, such 
as non-uniqueness of solutions (e.g. for L°° weak solutions of the MGo equation ll30l ). or norm explosion, 
i.e. arbitrary small data leads to an arbitrarily large solution in arbitrary small time, preventing the solution 
map to be continuous (cf. OH for the fluid equations, or PTll for dispersive equations). See Tao's book OTTl 
for a more detailed discussion. 

In this paper we prove that the solution map associated to the Cauchy problem for the MGo equation 
is not Lipschitz continuous with respect to perturbations in the initial data around a specific steady profile 
Oo, in the topology of a certain Sobolev space X. Hence the Cauchy problem is ill-posed in the sense of 
Hadamard. In order to prove this we first consider the instability of the MGo equation, when linearized 
about a particular steady state. We employ techniques from continued fractions in order to construct an 
unstable eigenvalue for the linearized operator. The use of continued fractions in a fluid stability problem 
was introduced by Meshalkin and Sinai [25 ] for the Navier-Stokes equations and later adapted for the Euler 
equations by Friedlander, Strauss, and Vishik lfl4ll . In contrast with the Navier-Stokes and Euler equations, 
where the linear operator generates a bounded semigroup, for the non-diffusive MG equation we construct 
eigenvalues which have arbitrarily large real part. This shows, in particular, that the linearized MGo operator 
does not generate a semigroup over Sobolev spaces. Ill-posedness in Sobolev spaces for the full nonlinear 
problem, then follows from perturbation arguments, which consist of showing that we cannot have a solution 
map that is Lipschitz with respect to the initial data. 

We emphasize that the mechanism giving the ill-posedness of the nonlinear MGo equations is not merely 
the order 1 derivative loss in the map h > U. Rather it is a combination of this derivative loss, with the 
anisotropy of the Fourier symbol of M in ( |1.5| >-( [L7| ), and the fact that the symbol of M is even. We note 
that the even nature of the symbol of M plays a central role in the proof of non-uniqueness for L°°-weak 
solutions to MGo given by Shvydkoy in [30], via methods from convex integration. In contrast, an example 
of an active scalar equation where the map 6 i->- U is unbounded, but given by an odd Fourier multiplier, 
is the limiting case of the inviscid modified SQG equation, introduced by Ohkitani ||2"8l . This equation was 
recently shown in Q to give a locally well-posed problem in Sobolev spaces (see also Remark |2~4| below). 

Lastly, we study the critically diffusive system ( |1.1| )-( [T3] ). We consider the evolution equation ( |1.1| > 
linearized about an arbitrary smooth steady state Go, and prove that if the associated linearized operator 
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has an unstable eigenvalue, then the full nonlinear AMG K system is Lyapunov unstable with respect to 
perturbations in the L 2 norm. The proof uses a bootstrap argument, a variant of which has been recently 
employed in a number of fluid contexts (see for example lfT2l[T3l ). Here the main difficulty arises due to the 
singular nature of U given by ( |1.3| >, which makes it challenging to control the growth of the nonlinearity. 
To overcome this obstacle we give a new a priori bound on a sub-critical Sobolev norm of 0, which is 
uniform in time. This bound does not follow directly from our earlier work [16], and we give the necessary 
details in the Appendix [A] We then turn to address the linear instability issue for MG K . Using a continued 
fraction construction, analogous to the treatment given in Section [2] for the case k = 0, we demonstrate the 
existence of unstable, but bounded eigenvalues. Combining this example with the result that linear implies 
nonlinear instability proven earlier for AMG K , we conclude that there exists a steady state around which 
the Moffatt model is nonlinearly unstable. Such an example produces strong growth in the magnetic field 
b(x, t), proportional to exp(t/«), which is consistent with the dynamo scenario. 

2. The non-diffusive equations 

In this section we consider the non-diffusive problem. First we address the AMGo equation, which we 
recall below 

d t e + U VG = (2.1) 
Uj = diTijG, V • U = 0. (2.2) 



2.1. Local well-posedness in analytic spaces for the AMGo equation. For a large class of equations 
which arise as singular limits in fluid mechanics (e.g. the Prandtl boundary layer equations, the hydrostatic 
Euler equations), the strength and structure of the nonlinearity prevents a local existence theory in spaces 
with finite order of smoothness, such as Sobolev spaces. However, if the derivative loss in the nonlinearity 
U • V© is of order at most one, both in U and in VO, it is possible to obtain the local existence and 
uniqueness of solutions in spaces of real-analytic functions, in the spirit of a Cauchy-Kowalewskaya result 
(cf. (H[T7l|22]| and references therein). In this section we prove that the AMGo equation falls in the category 
of singular equations for which we have analytic well-posedness, and we give a condition that prevents the 
breakdown of real-analyticity of the solution. Our first result in this direction is the following theorem. 

Theorem 2.1 (Analytic local well-posedness). Let 0(-,O) be real-analytic, with radius of analyticity at 
least To, and analytic norm Kq. Then, there exists T = T(tq, Kq) > and a unique real-analytic solution 
on [0, T) to the initial value problem associated to ( |2.1| )-( |2~l2"| ). 

Classically, a real-analytic function f(x) is characterized by a bound of the type \d a f(x)\ < K\a\W~\ a \, 
for all x and all multi-indices a, where K, r > are constants. The supremum over all values of r for which 
such a bound is available is the analyticity radius of /. In the periodic setting of this paper there exists an 
alternate characterization of real-analytic functions: / is real-analytic with radius proportional to r > if 
and only if the Fourier coefficients f(k) decay like exp(— r\k\) as \k\ — > oo. 



Proof of Theorem 2.1 As in lfTTll23ll24ll . we use the Gevrey-class real-analytic norm 



\\Q\\ 2 T = ||(-Ar/V(- A ) 1/2 9||i 2 = Y, \k\ 2r e 2T ^\Q(k)\ 2 , 

where r > denotes the analyticity radius, r > d/2 + 3/2 is the Sobolev exponent, and since we work in 
the mean-free setting we let = Z d \ {0}. We have the a priori estimate 

2|l|0|l' = r||(-A) 1 /4 ||2 + {u . v @ 5 (_A)-e 2 ^- A ) 1/2 e) = t\\(-A)^@\\ 2 T + K, (2.3) 



ILLAVELL-POSEDNESS RESULTS FOR THE MG SYSTEM 



5 



where (-, •) is the standard I? inner product. Using Plancherel's theorem, the nonlinear term 1Z may be 
written as 

n = i{2n) d Y u(j)-k@(k)\i\ 2r e 2T ^e(-i). 



Since the the matrix Ty consists of Fourier multipliers with bounded symbols, we obtain from ( |2.2[ ) that 
\U(j)\ < |j||9(j)| for all j G Z«f, and hence 

n < c Y lill fc l(lil r + |fcr)|e(j)|e T|, ' l |e(fc)|e r i*i|ine(i)| e T i , i 

<C Y CI J l"'"^ 172 1^| 3/2 -I" l^r - '" 172 | 3/2 ) I QCJ) le^l^'l |0C^)|e^l fc l |Zr H " :L / 2 |@(Z) |e^l z > (2.4) 

j+k=l;j,k,l&Zi 

for some positive constant C = C(r). In the above estimate we have used the triangle inequality \l\ < 
\j\ + |fc|, and |fc| 1/2 < \j\ 1/2 + |Z| 1/2 < 2\j\ 1/2 \l\ 1/2 since \ \l\ > 1. Since the right side of estimate ( gg) > 



is symmetric with respect to |j| and |fc|, we need to only give the bound for the case \j\ < \k\. Summing 
first in k, and using the Cauchy-Schwartz inequality we obtain from (|2.4|) that 



U < CU-A^QWl Y |j| 3/2 |@(j)|e T|il < C r ||(-A) 1/4 G|| 2 ||e|| r (2.5) 
since r > d/2 + 3/2, for some positive constant C r depending only on r and universal constants. From 



( |2.3| > and \2.5\ we obtain the a priori bound 

^ll (-'*)llr (t) < (r(i) + c r ||e(-,t)|| r(t) ) iK-A^eM)]! 2 ^. (2.6) 

Having arrived at ( |2.6| ), the argument is the same as in ||23ll24l . More precisely, we may let r(i) be decreas- 
ing and satisfy the ordinary differential equation 

t + 2C r K = (2.7) 

with initial condition r(0) = To, so that the right side of ( |2.6[ ) is negative, and therefore 

||e(-,t)|| TW <K = ||e(-,o)|| TO 

as long as r(t) > 0. These arguments prove the existence of a real-analytic solution 8(i) on [0, T*), where 
the maximal time of existence of the real-analytic solution is the time it takes the analyticity radius r(t) to 
reach 0, i.e. T* = tq/ {2C r Ko). This proof may be made formal using a standard Picard iteration argument. 
We omit further details. □ 



Whereas Theorem 2.1 guarantees the local in time existence and uniqueness of a real-analytic solution to 



( |2. 1| > — ( |272] >, it gives almost no intuition on what happens to the real-analytic solution at the time when the 



analyticity radius becomes 0. In Theorem 2.2 below we prove that if at time T a certain Sobolev norm of the 
solution is suitably controlled in terms of the initial data, then the real-analytic solution may be continued 
past T; therefore the breakdown of real-analyticity is controlled only by the initial data and only a Sobolev 
norm of the solution. 

Theorem 2.2 (Criterion for the breakdown of analyticity). Assume the initial data &(■, 0) is real-analytic 
with analyticity radius at least tq > 0, and analytic norm Kq > 0. For t > define the Fourier-Sobolev 
norm 

a{t) = || ((-A) 3 / 4 6(., S )) U = Y W 3/2 |e(M)|, (2.8) 
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and let 



A{t) 



a(s)ds. 



There exists a positive constant C r depending only onr > d/2 + 5/2, such that if at time T we have 

> A(T)exv(c r K exp(- A(t)) dtj (2.9) 



then there exists a real-analytic solution of the initial value problem associated to the AMGq equations on 
[0,T + 5], for some 5 > 0. 



The above analytic persistency criterion ( |2.9[ ) should be viewed in the context of numerical simulations 
of the AMGo equations: given Kq and tq, Theorem |2 . 2 | pro vides a numerically easy-to-track quantity which 
signals the time when the simulation becomes un-reliable/under-resolved. In particular, for most numeric 
simulations the initial data is entire (a finite trigonometric polynomial), and hence tq may be taken arbitrarily 



large. In this case condition (2.9 ) shows that if A(t) increases at a fast rate near T, the exponential factor on 



the right side of (Z9I becomes negligible compared to A(T), and so the breakdown of analyticity is due to 
the accumulation near T of the Sobolev norm. 



Proof of Theorem 2.2 With the notation of the proof of Theorem , we recall the estimate ( 2.3 1, i.e. 



1 d 

2dt 



101 



(2.10) 



The proofs differ in the way we estimate nonlinear term 1Z. From ( |2.4[ ), which is symmetric with respect to 
\j\ and |fc|, and using the bound e x < 1 + xe x which holds for all x > 0, 1Z may be bounded as 

n< c |j| 3/2 |e(j)|e r|j ' l |fcr +1/2 |e(fc)|e T|fe| |/r +1/2 |e(Z)|e T|i| 
j+k=l;j,k,leZi;\j\<\k\ 



< C r ||(-A) 1/4 e|| 2 ( J2\3\ 3/2 \&(j)\+r\\@\\ T 



(2.11) 



ifr > d/2 + 5/2, where C r > is a large enough constant depending only on r. Combining ( 2. 10 > with 
( 2. 1 1 [ >, and recalling the notation (2^8 ), we obtain the a priori estimate 



ld_ 
2dt 



\Q\\ 2 T <{t + C r a(t) + C r r||G|| r ) Wi-A^ef^. 



Therefore, if we choose r(t) decreasing fast enough so that 

f + 3C r a(t) + 2C r T\\e\\ T < 0, 

then from ( [2T2) ), using that ||6|| T < ||(-A) 1 / 4 e|| T , we obtain 

d 



dt 



\e\\ T + c r a(t)\\e\\ T < o 



(2.12) 



(2.13) 



(2.14) 



as long as r > and p,12| ) holds. Recalling that K = ||6(-, 0)|| T[) , it follows from ( |2.14D that ( |2.13| ) holds 
if we let r(t) solve the initial value problem 



f(t) + 3C r a(t) + 2C r K r(t) exp(-C r A(t)) = 0, 



(2.15) 
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with initial data r(0). This ordinary differential equation may be solved explicitly, yielding that 

T (t) = T exp (-2C r K J exp(-C r A(s))(is^ 

- 3C r J a{s) exp (-2C r K exp(-C r A(z)dz\ da. (2.16) 



After a short calculation, ( |2.16| > shows that that as long as 

ft 



3Ct 



> A(t) exp ^-2C r K J exp(-C r A(s))ds 



we have r(t) > 0, and therefore the real-analytic solution may be continued a bit past t, thereby concluding 
the proof of the theorem. □ 

Remark 2.3 (A proof in real variables). While for the AMGo equations set on T d the Gevrey-class norms 
introduced in ifTTTl are very convenient to work with, the proofs of both Theorem |2.1| and Theorem 2.2 may 
alternatively be given in the real variables. For this purpose, one may use the real-analytic norms introduced 
in E2l for the study of the hydrostatic Euler equations, i.e. one may work with 

linn ^ll^mi TH(|a| + l) 3 V-imaoi, rH-^lal + l) 3 



Id! " " T ^ " rial - 1)! 

a>0 11 a>i u 1 ; 



for some r > 0, where a G Ffjj is a multi-index. Closely following the arguments in ll22l . one may obtain 

-^ll©IU T < (f + C (l + r- 5 / 2 )||e||x T ) ||6||y T (2.17) 

for some constant Co > 0, as long as r > 0. The estimate ( 2.17| ) corresponds to the bounds (2^6) or ( 2.12 i 
from the periodic setting. A suitable choice of decreasing r then completes the proof in the real variables, as 
in the T d case. Although a real-variable proof is more natural, we have chosen here to present Fourier-based 
proofs of all results in this section for self-consistency of the paper, and transparency of the proofs. 



Remark 2.4 (Additional structure of the X^- implies local well-posedness in Sobolev spaces). If besides 

(T ij f,g) = (T ij g,f) (2.18) 



the fact that didjTijf = for all smooth /, the matrix {Tij}f =1 additionally satisfies 



for all smooth / and g, the AMGo equations are locally well-posed in Sobolev spaces. Note that if ( |2.18| ) 
holds, the operator i->- U is anti-symmetric (given by an odd Fourier symbol). In this case, the AMGo 
equations are locally well-posed in Sobolev spaces H s , with s large enough. The proof closely follows from 
the arguments of Chae, Constantin, Cordoba, Gancedo, and Wu 0. In particular, when estimating the H s 
norm of the most "dangerous" terms for closing the energy estimate are 

Tlx = (U- V(-A) s / 2 0,(-A) s / 2 0) and K 2 = {{-k) s/2 U ■ V6, (-A) s / 2 6). 

The first term 1Z\ is identically since U is divergence-free. Here one uses that didjTy = 0. For the second 
term one may write, using ( |2.18[ ) and integration by parts, 

n 2 = (diT^-Ay^e a,-e, (-A) s / 2 e) = -((-A) s / 2 e,^(a,e (-A) s / 2 e)) 

= -n 2 + ((-A) s / 2 0, [dj@, diTij}(-A) s / 2 @), (2.19) 



where the bracket [•, •] denotes a commutator. The second term on the right side of ( 2.19 1 is a Coifman- 
Meyer commutator, which is hence bounded as ||dj3j0||i,°o ||0||^. s . We omit further details and refer the 
reader to (21, where it was noted that the anti-symmetry of the map t- > U and the divergence-free nature 
if U may be used to close Sobolev estimates, via commutators. 
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2.2. Hadamard ill-posedness the linearized MGo equations. The classical approach to Hadamard ill- 
posedness for partial differential equations arising in fluid dynamics (see for instance lfT7l[T8l[T9ll2"0l|2"9l ) is 
to first linearize the equations about a suitable steady state Qq, such that the linear operator L obtained has 
eigenvalues with arbitrarily large real part, in the unstable region. In particular, this very strong instability 
shows that the linear operator L does not generate a semi-group: 

Theorem 2.5 (Linear ill-posedness). The Cauchy problem associated to the linear evolution 

d t 9 = L9 (2.20) 

where the linear operator L and the steady state Go are given by 

L9(x, t) = -M 3 9(x, t) d 3 @ (x 3 ) (2.21) 

©0(^3) = asin(mx3) (2.22) 

for some a > 0, and integer m > 0, is ill-posed in the sense of Hadamard over L 2 . More precisely, for any 
T > and any K > 0, there exists a real-analytic initial data 9(0) such that the Cauchy problem associated 
to (2.20 1— ( 2.22 1 has no solution 9 £ L°°(0, T; L 2 ) satisfying 



sup \\6(;t)\\ L 2 <K\\9(0)\\ Y 
te(o,T) 



(2.23) 



where Y is any Sobolev space embedded in L 2 . 

In order to prove Theorem|2.5| we need a more detailed spectral analysis of the linear operator L defined 



in ( 2.21 1. We consider the simplest possible steady state, namely Uq = 0, ©0 = ^(^3), for some T-periodic 
function F, with F(x 3 ) dx 3 = 0. We denote by 9 the perturbation — ©o- Linearizing ( |2.1[ ) about this 
steady state gives 



d t e + u 3 F'( X3 ) = o, 

where u 3 = M38 is given explicitly by ( |1.7| > via its Fourier coefficients 



u 3 (k) 



An 2 k 2 \k\ 2 + n 2 kl 



6{k) 



(2.24) 



(2.25) 



where we have denoted the physical parameter \i = /3 2 / 77 > 0. It is then natural to consider F to be a 
27r -periodic function consisting just of one harmonic, namely 

F(x 3 ) = asin(m3;3) (2.26) 

for a fixed amplitude a > 0, and an integer modulation m > 0. Note that F'(x 3 ) = mocos(ms3). 

Since the linear evolution ( 2.24 1 does not act on the x\ and X2 variables, it is natural to pick initial data 
#(-, 0) consisting of only one Fourier mode in x\ and x-i- We hence look for solutions to ( |2.24[ ) of the form 



6(x, t) = e at sin(/cixi) sin(/c2^2) Cn sm ( nx 3) 



(2.27) 



n>l 



for some fixed integers k\ and fo- In ( 2.27 1, we have chosen a sine series for convenience, since f Q * 9 dx 
0. Therefore, by ( |2.25 1 we also have 



u 3 (x, t) = e at sin(&i2;i) sin^^) c ' r ' 



nk 2 (kf + 



n>l 



4n 2 n 2 (k 2 + k 2 + n 2 ) + fi 2 k' : 



sin(nx3). (2.28) 



After inserting ( 2.27 1 and ( 2.28 1 into ( 2.24[ ), and dividing by e at sin(fcixi) sin(A;2a;2), we obtain 



a c n sin(nx3) + am cos(mx 3 ) c n 



p,k 2 (kf + k 2 ) 



n>l 



n>l 



An 2 n 2 (k 2 + k 2 + n 2 ) + p 2 ^ 



sin(nx3) = 0. 



(2.29) 
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To simplify the interaction between cos(mxs) and sin(nxs), for the value of m fixed in ( 2.26 ), we make the 
additional ansatz 



c n = 0, whenever n is not an integer multiple of m, 
c n = Cmp =-c p , whenever n = mp, for some integer p > 1. 
Inserting ( [230l >- p3T] ) above into ( [2T29] ) gives 

a ^2 Cp sin(mp X3) + — ( sin(m(p + 1) X3) + sm(m(p — 1) X3) 

p>i p>i p 

where for ease of notation, since k\ and k2 are fixed throughout, we have denoted 



ft. 



2 3 Q 2 {mp) 2 (k 2 + k 2 2 + {mp) 2 ) + 2^i 2 k\ 
a\xmk\{k\ + k 2 ) 



(2.30) 
(2.31) 



(2.32) 



(2.33) 



for any p > 1, where a, m > 1 are fixed as in ( |2.26| ). The essential feature of the a p coefficients is that they 
grow as p 4 when p — > 00. The above equation ( |2.32[ ) gives the recurrence relation for the sequence c p , in 
terms of the given sequence a p 



(Tc p H 1 

C-2 



0, for all p > 2, 



CJCl + 



«2 



0, for p = 1. 



Letting r] p = (c p a p -i) / '(c p -ia p ) , we obtain 



+ Vp+i + 



?7p 



0, p > 2 



and therefore for any p > 2 we have 



ffai +772 = 0, p = 1, 



-1 



-1 



(2.34) 
(2.35) 

(2.36) 
(2.37) 

(2.38) 



P iP^ x P cra p+l +r) p+2 

and we obtain rj p as a continued fraction. Equating with the case p = 2 gives an equation in a, namely 

. (2.39) 



1 



aa2 



aa :j - 



The goal is to find a real, positive, solution a of (239), since then the solution 8(t) grows exponentially in 
time like exp(crt) (cf. ( 2.27j >). The following lemma states the existence of such a solution. 



Lemma 2.6. Let a p be defined for all p > 1 as in (2.33 ), where the positive integers m, a, k\, ^2 are fixed, 
and n,Q > are fixed physical parameters. Then, there exists a real solution <r* > to ( 2.39| ), and 
moreover we have a* > 1/ - v /ai02. 



Proof of Lemma 2.6 For ease of notation we define the functions F p and G p as 

W = 



*p+2- 



G p (a) 



a 2 a 2 - 4 



aa p + J<j 2 a 2 — 4 



(2.40) 



(2.41) 
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for all p > 2, and all a > do = 2/«2- -Pp is well-defined and smooth, except for a set of points on the real 
axis with a < gq. For the rest of the proof of this lemma we will always assume that a > o"o > 2 / a p for 
any p > 2. Note that G v satisfies 



G p (a) 



1 



aa p - G p (a) 



oar 



for all for p > 2. Since the sequence a p is strictly increasing in p, we have that 

G 2 ((t) > G 3 (a) > Gi(a)... > 0. (2.42) 
Since linip^oo a p = oo, we have that linip^oo G p (a) = for every fixed a. Lastly, for all p > 2 we have 



oa p > G p+ \(a 



(2.43) 



for all a > cr . From ( |2.41[ ) and ( |2.42[ ) we obtain that cra2 — ^(cr) > cjq2 — G2(er) > and therefore 

1 1 



G 2 (<r) 



> 



cra2 — G2(a) cra-2 — G3(cr) 
Similarly, to the above inequality, using ( |2.42[ ) and ( |2.43| >, we obtain 

1 



> 0. 



< aa 2 — Gs(a) = aa 2 



aa 3 - G 3 (a) 



< aa.2 



so that 



G 2 {o) > 



1 



> 



aa<2 — Gs(cr) cra 2 



cr«3 — G±{cr) ' 



> 



<T03— Gi(cr) 



holds for all a > do- An inductive argument then gives that 

„ , , 1 1 



aa>2 — Gs(a) 0012 V< 7 \ 

y ' z 0-03— GaKP) 



oa.2 



F 2 (a) > (2.44) 



0-013—.. 



for all p > 2, and a.e. a > ctq. Repeating the above construction argument we obtain that < ^3(0") < 
Gs(a), and hence by ( |2.43[ ) we get < aa.2 — Gs(a) < aa2 — F 3 (a) < aa2, so that 



F 2 (<r) 



> 



aa2 — F 3 (a) aa2 
Therefore, it follows that there exists a real cr* > ao = 2/ 02 such that 

a*ai = F 2 (a*), 

and moreover, from ( |2.44[ ) and ( |2.45| > we obtain the estimate (see Figure[T]below) 



1 



which concludes the proof of the lemma. 



< cr* < 



\J a\a2 — a\ 



(2.45) 



(2.46) 



(2.47) 



□ 



From Lemma 2.6 we therefore obtain a solution cr* to ( |2.39[ ). Since the coefficients a p are given 
(cf. ( 2.33 1), the recursion relations ( 2.36 l-( 2.37[ ) uniquely defines the values of rj p . We choose 



ci = a\, 

c p = a p r} v r) p _x ■■■V2, for all p > 2. 



(2.48) 
(2.49) 
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FIGURE 1. Diagram showing how estimates (2.44 1 and (2.45 1 give the bound (2.47 1 



This sequence satisfies the recursion relations ( |2.34| )- p.35| ) by construction. We claim that the above defined 
coefficients c p decay very fast as p — > oo. To see this, from ( |2.38[ ) we observe that 

-1 -1 _ -1 

rjp — ■ — | — j — ... — —rp^a*). 



a*a p + r]p + i a*a % 



l 

cr t a p+ i+r] p+ 2 



v 



< T * a v +2+ r lp+3 



Moreover, similarly to ( |2.44| )-( |2.45[ ) one may prove that G p {a) > F p (a) > l/(aa p ) for all a > oq, and 
hence 

-2 -1 

< Vp < — . (2.50) 



cr*a p + y a%a l - 4 



cr*a 



v 



Moreover, from ( |2.33| > we have a p 

0(l/p 4 ) as p -> oo. Using (|230]>, it follows from (gggb that 



0(p ) as p — >• oo, and we thus from ( |2.50[ ) above we obtain r/ p 



lim c = 



and this convergence is very fast, namely of the order of C p /((p — l)!) 4 as p — > oo, for some positive 
constant C = C(cj*, /j,, a, m, k±, kz). Therefore, the solution 9(x, t) of ( |2.26 1, defined by ( |2.27 1 lies in 
any Sobolev space, it is C°° smooth, and is even real-analytic. In summary, we have proven: 

Lemma 2.7. Fix integers a, m > 1 arcc? a Sobolev space Y C L 2 . For any integer j > m there exists a 
C°° smooth initial datum 9^\-,0), with 0)||y = 1, and a C°° smooth function 9^\x, t) solving the 

initial value problem associated to the linearized MGo equations ( 2.20[ )-( [2.22[ ), such that 

) (2.51) 

for all t > 0, where C a)TO)jU) n > g/verc ( |2.52| ) below. 

Proof of Lemma 2.7 Recalling the definition of a p (cf. ( |2.33| >), we have that the function 9(-, t) given by 
( |2.27| ), with coefficients c p given by ( |2.48| )- p.49[ ), grows at the rate of exp(cr*t), where is such that 

a\imk\(k\ + k\) 



2 5 Q 2 m 2 (kl + k 2 2 + Am 2 ) + 2^ 2 fc| ' 
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Letting (k\, k 2 ) £ Z 2 be such that &| = k% = j > m > 1, we obtain that for any a, m, and j with j > m, 
there exists a solution <r* of (239 1 bounded from below as 

a\im 



cr*> J 



2 8 n 2 m 2 + 2p? 



jCa. 



concluding the proof of the lemma. 



(2.52) 
□ 



In order to prove the severe type of ill-posedness stated in Theorem 2.5 we need to show that it is not 
the uniqueness of solutions to the Cauchy problem associated to ( |2.20 i which fails. Indeed, we have the 
following result regarding the uniqueness of smooth solutions to the linear equation. 

Proposition 2.8 (Uniqueness for the linearized equations). Let 9 £ L°°(0,T;L 2 ) be a solution of the 
Cauchy problem associated to ( 2.20 1— ( 2.22 ), with initial data 9(-, 0) = 0. Then for any t £ (0, T) we have 

9(;t) = 0. 

Proof of Proposition 2.8 Since 9 £ Lf°L 2 , it can be written as the sum of its Fourier series 

0{x,t) = Y 0(k,t)e lk - x . 
fcez 3 



Taking the Fourier transform of ( |2.20| ), and noting that c^Qo is a function just of x 3 , we obtain (suppressing 
the t dependence) 



d t 9(k) + ]T 



AWn 2 {kl + k 2 2 +n 2 ) + ^k* 



i, k 2 , n)i(k 3 - n)e (fc 3 - n) 







nez\{o,fc 3 } 

where we denoted the Fourier series coefficients of Bo by Qq(ti), for all n £ Z \ {0}. Multiplying the above 
equation by 6(k) and summing over k% £ Z \ {0}, we obtain the bound 



1 d 

2dt 



\0(h,k 2 ,k 3 )f 



fc 3 ez\{o} 



X/ 40 2 n 2 
fc 3 ez\{o} nez\{o,fc 3 } 



< 



< 



lik 



' "9 3 6o||l2 



40 2 

flk 2 7T 2 




I ^(fci , k 2 , n) 1 1 k 3 - n\ 1 9 ( k 3 - n) \ \ 9(h , k 2 , k 3 ) \ 

1/2 



\0(h,k 2 ,h)f 



E 



|0(A;i,fc2,n)| 



1 



fie 



" i>v '' fc 3 ez\{o} 



(2.53) 



where in the second inequality we reversed the order of summation and used the Cauchy-Schwartz inequality 
in the k% summation, and in the third inequality we used that 



n>l 



n 



■ vr 4 /90. Note that in ( [233] ) we do 
not sum in k 2 , nor in k\. Using Gronwall's inequality, and the initial condition 9(0, k) = for all k £ Z 3 , 

(2.54) 



we obtain from (12.531) 



Y \0(t,k 1: k 2 ,k 3 
fc 3 ez\{o} 







for any fixed ki,k 2 £ Z, and all t £ (0,T). Since k\ and ^2 are arbitrary, it follows from ( |2.54[ ) that 
9(t, k) = for all k £ Z 3 , with k 3 ^ 0, and all t £ (0, T), concluding the proof of the proposition. □ 



We now have all the ingredients needed to prove Theorem 2.5 The guiding principle behind this result is 
that the existence of eigenfunctions for L, with arbitrarily large eigenvalues in the unstable region, prevents 
the existence of a continuous solution map, and in particular the problem is hence ill-posed in the sense of 
Hadamard. 
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Proof of Theorem [23] Let T > and K > be arbitrary, and let #(-, 0) 6 y C I 2 . Assume by contradic- 
tion that there exists a solution 9 £ L°°(0, T; L 2 ) of the Cauchy problem associated to ( 2.20| )-( |2l22 1 such 
that (12.231) holds, i.e. 



sup ||0(-,t)|| L 2 <K 
te[o,T] 



(2.55) 



Note that this solution also satisfies dtO £ L°° (0, T; ) so that 9 is weakly continuous with values in L 



On the other hand, for any j > m, from Lemma 



2.7 



we obtain the existence of a smooth solution 9^' of 



( |2T20l ), with ||0(-?)(-,O)||y = 1, which by ((23T) satisfies 

||^(-,i)|| L2 >exp(iTC Wj n/2), 



(2.56) 



for all t £ (T/2, T). Moreover, this solution is unique in LfL 2 . by Proposition 2.8 Letting j be a suffi- 
ciently large integer, for instance choosing 



J 



21og(l + iT) 
max < — - — — , m 



TC 'a,m,fj,,Q 

we obtain from ( |2.56| > a contradiction with ( |2.55| ), thereby proving the theorem. 



□ 



2.3. Hadamard ill-posedness in Sobolev spaces for the nonlinear MGo equation. We first make precise 
the definition of Lipschitz well-posedness for the full nonlinear MGo equation (see |[20l Definition 1.1]). 



Definition 2.9 (Lipschitz local well-posedness). Let Y C X c 

problem for the MGq equation 

d t e + u- ve 



W ' be Banach spaces. The Cauchy 

(2.57) 
(2.58) 

is locally Lipschitz. (X, Y) well-posed, if there exist continuous functions T,K : [0, oo) 2 — > (0, oo), the 
time of existence and the Lipschitz constant, so that for every pair of initial data Q^\-, 0), Q^ 2 \-, 0) £ Y 
there exist unique solutions £ L°°(0, T; X) of the initial value problem associated to the MGq 



V-U = 0, U. 



j 







equation ( 2.57 )-( 2.58 1, that additionally satisfy 

\\Q( 1 )(.,t)-e^(;t)\\ X <K\\e^\;0)-Q^\;0)\\ Y (2.59) 

foreveryte [0,T], T = T(||G( 1 )(-, 0)||y , ||0< 2 >(-, 0)||y) am/if = iT(||e( 1 )(-, 0)||y , ||G( 2 )(-, 0) ||y). 

If G( 2 )(-, t) = and X = Y, the above definition recovers the usual definition of local well-posedness 



with a continuous solution map. However, Defintion 2.9 allows the solution map to lose regularity. In EOl 
the above defined well-posedness is referred to as weak well-posedness. It has been pointed out to us by 
Benjamin Texier that a loss in regularity is usually needed in order to obtain Lipschitz continuity of the 
solution map for quasi-linear equations, and hence the typical spaces (X, Y) we have in mind are X = H s , 
and Y = H s+1 , with s > 1 + d/4. For the purpose of our ill-posedness result, we shall let Q^ 2 \x, t) be 
the steady state 8(^3) introduced earlier in ( |2.22| ). We consider X to be a Sobolev space with high enough 
regularity so that <9 t £ L°°(0, T; L 2 ), which implies that is weakly continuous on [0, T] with values in 
X, making sense of the initial value problem associated to ( 2.57[ )- p.58| ). 
The main result of this subsection is the following theorem. 

Theorem 2.10 (Nonlinear ill-posedness in Sobolev spaces). The MGq equations ar e loc ally Lipschitz 
(X, Y) ill-posed in Sobolev spaces Y C X embedded in W 1,4 , in the sense of Definition 



2.9 



above. 



The proof of Theorem 2.10 follows from the strong linear ill-posedness obtained in Theorem 2.5 



com- 



bined with the uniqueness of solutions to the linearized equations proven in Proposition 2.8 and a fairly 
generic perturbative argument (cf. OT1 pp. 183]). This program has been successfully used in the context of 
the hydrostatic Euler and Navier-Stokes equations ll29l . or the Prandtl equations UHdHHO]], an( ^ others. 
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Proof of Theorem 2.10 Let Y C X C W 1,4 be as in the statement of the theorem. Since X embeds in 
H , the linearized operator L@ = — M36 c^Go maps X continuously into L 2 , and since X C W 1,4 , the 
nonlinearity iVG = -MjQ djQ be bounded as ||./V0|| L 2 < ||V6||| 4 < C\\ 9 ||^-, for some constant C > 0. 
For instance, one may consider X = H s and Y = H s+1 , where s > 7/4 in three dimensions. 

Fix the steady state 60(^3) G Y, as given by ( |2.22| ). Also, fix a smooth function iJjq G Y, normal- 
ized to have HV'oIIy = 1, to be chosen precisely later. The proof is by contradiction. Assume that the 
Cauchy problem for the MGo equation ( |2.57| )-( |2.58[ ) were Lipschitz locally well-posed in (X, Y). Consider 
9( 2 )(:z,0) = 60(2:3), so that e^ 2 \x,t) = o (x 3 ) for any t > 0. Also let 

for every < e < ||0o||y. To simplify notation we write e instead of Q^ 1 ' 6 ). By definition 2.9 for 
every e as before there exists a positive time T = T(||Bo||y, ||@ e ||y) and a positive Lipschitz constant 
K = (||Oo||y, ||© e ||y) such that by ( |2.59| ) and the choice of ipo we have 

||e e (-,t) - e (-)IU < (2.60) 

foralUG [0,T]. We note that since ||B e (-,0)||y < ||6 ||y + e < 2||6 ||y, due to the continuity of T and 
K with respect to the second coordinate, we may choose K = if(||0o||y) > and T = T(||0o||y) > 
independent of e G (0, ||0 o ||y), such that ( |2.60| > holds on [0, T}. 

The main idea is to write the solution e as a perturbation of 0o, i.e. define 

e e (x,t) = e e (x,t) -e (x 3 ) 

for all t G [0, T] and all e as before. The initial value problem solved by 9 e is therefore 

d t e e = L9 e + Nd e 

where we recall that the linear and nonlinear operators are given by 

L0 e = -M 3 9 e d 3 e 



N9 e 



-Ma 



and Mj is defined by its Fourier symbol cf. ( |1.5[ )-( [l~7j ) for j G {1, 2, 3}. It follows from ( 2.60 ) that the 
function 

^(x,t) = 9 e (x,t)/e 

is uniformly bounded with respect to e in L°°(0, T; X). Therefore, there exists a function ip, the weak-* 
limit of ip e in L°°(0, T; X). Note that tp e solves the Cauchy problem 

d t ^ = W + eN(?p € ) (2.61) 

tl> e (;0)=ik- (2-62) 

Due to the choice of X, we have the bound 

\\N^\\ L 2<C\\r\\x<CK 2 , (2.63) 

and from ( |2.61[ ) we obtain that dtip e is uniformly bounded with respect to e in L°°(0, T; L 2 ). Therefore the 
convergence ip € ^ ip is strong when measured in L 2 . Sending e to in ( |2.61[ ), and using ( |2.63[ ), it follows 
that 



^j(x,0) 



Lip 

■00 



holds in L°°(0, T; L 2 ), and this solution is unique due to Proposition |2.8| In addition, the solution tp inherits 
from ( |2.60| ) the upper bound 

U(-,t)\\ L2 <K (2.64) 
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for all t G [0, T}. But this is a contradiction with Theorem 2.5 Indeed, due to the existence of eigenfunctions 
for the linearized operator with arbitrarily large eigenvalues, one may choose ipo (as in Lemma 2.7 ) to yield 
a large enough eigenvalue so that in time T/2 the solution grows to have L 2 norm larger than 2K, therefore 



contradicting ( 2.64 >. □ 



3. The critically diffusive equations 

In this section we consider the dissipative equations, where k > is fixed. First we prove for the 
AMG K equations that linear implies nonlinear instability, via a bootstrap argument that is a combination of 
arguments given in lTT2l[T3l . and a new a priori bound obtained in the appendix. We then turn to the specific 
MG K equations, for which we explicitly construct an unstable eigenvalue of the linear operator obtained 
when linearizing about a specific steady state. 

3.1. Linear instability implies nonlinear instability. For a smooth source term S G C°° we consider an 
arbitrary smooth steady state @o G C°° of ( |l.l| )-( fL"2"| ), i.e. a solution of 

U • V6o = S + kAOq, (3.1) 

where the divergence-free steady state velocity is given by 

U 0j = diTijQ . (3.2) 

Writing the AMG K equations < | 1 . 1 1 > — < fT72| > in perturbation form, we obtain that the evolution of the perturba- 
tion buoyancy 9 = — 0o is governed by the initial value problem associated to 

d t 9 = L9 + N9, (3.3) 

where the dissipative linear operator L is defined as 

L6 = -U -V9-U- V0 O + kA6 (3.4) 

with Uj = diTijO = Uj — Uoj, and the nonlinear operator N is given by 

NO = -u-V6 = -V ■ (u6). (3.5) 



In Section 3.2 below we prove that the linear operator L has eigenvalues in the unstable region in the specific 



case matTy = -5i(-A) -1 Mj. Here we recall a suitable version of Lyapunov stability (cf. llT2l[T3l[l4l ). 

Definition 3.1 (Nonlinear instability). Let (X, Z) be a pair ofBanach spaces. A steady state ©o is called 
(X, Z) nonlinearly stable if for any p > 0, there exists p > so that if0(-, 0) G X and \\9(-, 0)\\z < p, then 
we have 



(i) there exists a global in time solution 6 to the initial value problem (3.3 1— ( 3.5 1 with 9 G C( [0, oo) ; X); 

(ii) and we have the bound \\9(- ,t)\\z < pfora.e. t G [0, oo). 

An equilibrium ©o that is not stable ( in the above sense) is called Lyapunov unstable. 

The Banach space X is the space where a local existence theorem for the nonlinear equations is available, 
while Z is the space where the spectrum of the linear operator is analyzed, and where the instability is 
measured. For the AMG K equations we measure the instability in the natural space Z = L 2 (T d ), while for 
the space of local existence we may for instance consider X = H s (T d ), for any s > d/2. Note that H s 
with s > d/2 is subcritical for the natural scaling of the equations, and hence obtaining the local existence 
of a unique H s solution is not difficult. Now we are ready to formulate the main result of this section. 

Theorem 3.2 (Linear implies nonlinear instability). Let ©o be a smooth, mean zero steady state solution 
of the critically dissipative AMG K equations, i.e., it solves < |3- 1| >-( [3T2] )- If the associated linear operator L, as 
defined in ( |3.4[ ), has spectrum in the unstable region, then the steady state is (H s , L 2 ) Lyapunov nonlinearly 
unstable, for arbitrary s > d/2. 
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In order to prove Theorem |3.2| fix a linearly unstable smooth eigenfunction <j), with eigenvalue of maximal 
real part A. Fix a parameter 8 G (0, C\), where C\ > is to be determined. It is convenient to consider the 
operator 

L s = L - (A + 5)1 (3.6) 

which is a shift of L that moves the spectrum to the left of the line hnz = — 5 in the complex plane, so 
that the resolvent of Ls contains the full right-half of the complex plane. Similarly to the Navier-Stokes 
equations [12J, the operator L$ is a bounded lower order perturbation of the Laplacian (since both u and 
X79 only lose one derivative with respect to 6). Therefore, L$ generates a bounded (due to the shift) analytic 
semigroup e Lst over any Sobolev space, and in particular over L 2 . In consequence, we have 

WLle^U^ < Mf« (3.7) 

for any a > 0, and for any t > 0, where M is a positive constant. Also, since Ls is a bounded perturbation 
of the Laplacian, similarly to lPT2l for a bounded perturbation of the Stokes operator, we have 

II V"IU a ^ ^ll(-A)- a /|| L2 < C\\f\\ L2d/(d+4a) (3.8) 

for some positive constant C, a G [0, <i/4), and all smooth, mean-free functions /. Note that in the last 
inequality of the above estimate we have also used the classical Riesz potential estimate. 



We now turn to proving Theorem 3.2 In order to achieve this, it is sufficient to prove that the trivial 
solution 9 = of (3.3 1 is (H S ,L 2 ) nonlinearly unstable. This is achieved by considering a family of 



solutions 6 e to ( |3.3| ), with suitable size-e initial data, i.e. 

d t 6 e = L6 e + N(6 e ), (3.9) 
9 e \ t =o = e(f>, (3.10) 

where <fi is an eigenfunction of L associated with the eigenvalue of maximal positive real part A. We will 



prove the following proposition which clearly implies Theorem 3.2 



Proposition 3.3. There exist positive constants C* and e* < 1, such that for every e G (0, e*), there exists 
T t >bwith\\e € {T e )\\ L 2>C*. 

In summary, the above proposition states that an 0(e) perturbation of the trivial solution leads in finite 



time to an 0(1) perturbation, where e is arbitrarily small. We remark that if 6 e solves ( |3-9[ >— ( |3TT0] >, then 
9 e = G + 9 € solves (|T7T|)— (|TT2]> , with initial data 6 + ecp £ C°°. The main difficulty is to control the 
nonlinearity, so that the exponential growth arising from the linear part of the equations is not canceled by 
nonlinear effects. As in |[T3l for the critical SQG equations, this is achieved by making use of an a priori 
bound on a sub-critical norm of 0. In the case of the AMG K equations the following estimate turns out to 



be sufficient for proving Proposition 3.3 



Lemma 3.4. There exists a positive constant C* = C*(Qq, (f),S,€*) > such that 

||A0 e (-,i)|| i2 < C*, for all e G (0,e*) (3.11) 

and arbitrary t > 0. 

The proof of Lemma |3.4| does not follow directly from our earlier works lfl31 [T6l on the global well- 
posedness and higher regularity of the AMG K system, and we give the proof in Appendix |A|below. 



Proof of Proposition \3.3\ The proof is a hybrid of the proofs for the corresponding results for the Navier- 
Stokes lTT2l and the critically dissipative SQG equation lfT3l . It follows the bootstrap argument considered 
previously for the equations of fluid mechanics, see for example HI [12] [131 . The argument consists of 
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writing the mild formulation of (3.9 >— ( 3. 10 >, that is 

9 e (t) = ee Lt (j)+ [ e L ^N6 e (s)ds 
Jo 



ee Lt (b + 



A X+8) (t-s) e L s (t-s) NQ e 



ee Lt (b + 



,(X+6)(t-s) L a e L e (t-s) L ja Nd e^ dg (3 12 ) 

where a G (1/2, 1) is arbitrary, to be chosen later. Applying the L 2 norm to ( 3. 12 > we obtain from (3^1 and 
( [3~8] > (with d = 3), that 

\ L2 <C <t> ee xt + r e ( A + 5 )(*- s )||L^('- s )|| £(i2) ||L,- Q iVr( S )|| L2 ^ 



< C 6 ee M + CM 



o (X+S)(t-s)_ 



o 



\L7<*V-(u*(s)e*(s))\\ L 2ds 



< C 6 ee M + CM i c 
Ji 

rl 

< C 6 ee M + CM 



/' 

Jo 



(\+8)(t-s)_ 



(t 



\L 



1/2-Q 



u^s)9*(s))\\ L 2 ds 







(t - s. 

By the Holder inequality and interpolation of Lebesgue spaces, we have 

\\u e (s)6 e (s)\\ La/(ia+1) < C||0 e || i2 ||u e || i6 /(4c*-2) 



u e (s)6 e (s)\\ L 6/(4. a +i) ds. 



(3.13) 



< c\\e 



e na+l/4||^ e ||7/4-a 



(3.14) 



\ L 2 \\^ \\ L 2 

where the last inequality holds for any a G (3/4, 1), since in three dimensions L 6 is the largest Lebesgue 
space that embeds in H l . From here on we may fix a value of a G (3/4, 1) and for ease of notation we 
let henceforth a = 7/8. Therefore, inserting estimate ( 3.14 i into the bound < \3.13 1, and recalling ( |3.11| >, we 
obtain 

rt 



L 2 



< C^ee xt + 



C(C*) 7/8 [ 
Jo 



,(A+*)(t-«)_ 



1 



(t - s y/ s 



\% 8 



ds. 



(3.15) 



Since the exponent of the L 2 norm of 9 e in the above integral is larger than 1, one may conclude the 
proof of the lemma, by following mutatis mutandis the proof we gave for the critical SQG equation in lfl3l 
Proposition 4.1]. For the sake of completeness we give here a few details. For R > C^ := \\4>\\l 2 to be 
chosen later, let T = T(R, e) be the maximal time such that 



L 2<eRe xt , forte [0,T]. 



(3.16) 



Clearly T G (0, oo] due to the strong continuity in L 2 of t h- > 9{t) and the chosen initial condition. There- 
fore, letting 5 G (0, A/8), we obtain from ( |3.15| > and ( |3.16| > that on [0, T] we have 

< Qee At + Ci (eRe xt ) 9/8 



(3.17) 



where C\ = C\{C* , A, 5) > is a sufficiently large constant. Using the definition ( 3.16 ) and the above 
estimate, we obtain the following bound on T: 



18 



SUSAN FRIEDLANDER AND VLAD VICOL 



since we have chosen R > C$. Therefore we obtain that 

log(C 2 /e) 



T > T f 



A 



(3.19) 



To find a lower bound on ||6 )e (T e )|| i 2 one again uses the Duhamel formula bound ( 3.17 ) and the triangle 
inequality 



\e e (T e 



\L* 



XT. 



Ci eRe 



XT, 



9/8 



C0C2 — CiR^^C^ 8 — C2 (2C^ — R) —: C* 



thereby concluding the proof of Proposition 3.3 upon letting Ca, < R < 20s- 



□ 



3.2. Example of an unstable eigenvalue for the MG K equations. In this section we construct a linearly 
unstable eigenvalue for the critically dissipative MG K system. This shows that the assumptions of Theo- 
rem 



3.2 are satisfied for instance by the specific steady state constructed here. 



As in the non-diffusive case, we consider the steady state Go = F{x 3 ) and Uq = 0, where the source term 
is S = — kA0o. The example we work with is similar to Section[2j namely we let F(x^) = asm^mxs), 
with S = nam 2 sin(mx3), for some positive integers a,m > 1 to be determined below, in terms of k. We 
consider the linear evolution of the perturbation temperature 9 = — Go, and as in ( 2.27 ) we make the 
assumption that 6 is of the form 



9(x, t) = e at sin(&ixi) sin(/c2X2) / c n sin(nx3 



n>l 



which by ( |1.4[ ) implies that the perturbation velocity is 

uz(x, t) = e at sin(&ia?i) sin(A;22;2) A 



fik^kf + fc 2 ) 



40, 2 n 2 (k 2 + ki + n 2 
n>± 1 z 

Inserting ( 3.20[ )-( |3.21[ ) into the linearized MG K equation, we obtain 

+ n(k 2 + k 2 + n 2 ))c n sin(nx3) 



+ H 2 k 



—g sm(nx 3 ) 



(3.20) 



(3.21) 



n>l 



+ am cos (771x3) 



n>l 



l 4n 2 n 2 (k 2 + k 2 + n 2 ) + v 2 kl 



sm(nx3 1 



0. 



We now follow the construction of an unstable eigenvalue given in Section 2.2 for the case where k 
The difference in the analysis is that a is replaced by a p , where 

a + K(k\ + kl + m 2 p 2 ) 



0. 



(3.22) 



for all p > 1. A suitable modification of the method of continued fractions (see also E51 ) again produces 
a characteristic equation that gives a modification of the lower bound ( 2.47 ) on the root a* . Namely, in the 
diffusive case we have 

a\xmk\{k\ + k 2 ) 



a* > cj« > 



2 5 n 2 m 2 (kl + k\ + 4m 2 ) + 2y?k\ 



K(kf + k^ + 4m 



(3.23) 



Therefore, for at most finitely many values of k\ and k2 we have a* > 0, and we may choose a large enough 
in terms of k, m, [i and Q, so that there exists at least one pair {k\, /C2) giving rise to unstable eigenvalues. 

More precisely, we observe that when < /t < 1, which is the range relevant to the Earth's fluid core, 
the maximum value of the lower bound given in ( 3.23 1 as a function of the wave numbers k\ and k2 occurs 
when 
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for < k >C 1. Inserting this into (3.23 1, the lower bound on a* becomes 



•* > 2m -\ <3 ' 24) 

whenev er < k <C 1. Thus, the limit k — > is consistent with the non-diffusive problem studied in 
where it was shown that there exist eigenvalues a that grow like k\ on the curves k\ ^ k\, as 



2.2 



Section 

k\ — > oo. 

The classical geodynamo problem of magnetic field generation via the flow of an electrically conducting 
fluid is closely connected with the existence of unstable eigenvalues (see, for example, EoTl ). We have 
shown that the MG K equations are globally well-posed (cf. lfT5l ) and that the linearized equations permit 
solutions that grow exponentially in time. We note that this growth in time for 9(x, t) implies the growth of 



the perturbation magnetic field b(x, t) via the relation ( |1.10| >. The exponential growth rate is very rapid, i.e. 
of order 1/k, for very small buoyancy diffusivity k. 

Acknowledgements. The work of S.F. is supported in part by the NSF grant DMS 0803268. We are grateful 
to Benjamin Texier for pointing out that the definition of Lipschitz local well-posedness with X = Y would 
be too restrictive in the non-diffusive setting. 

Appendix A. On the higher regularity of the forced critical AMG k equations 



The goal of this section is to give a proof of Lemma 3.4 The main difficulty is to show that the constant 



C* in ( |3.1 1[ ) may be taken independently of time. Due to the triangle inequality 

||A0 e || L 2 < ||Ae e || L 2 + ||A6o|| L 2, (A.l) 



in order to prove Lemma 3.4 it is sufficient to consider bounds for the initial value problem associated to 

d t e e + u e ■ ve e - KAe e = s (A.2) 

V • U e = 0, Uj = diTijQ 6 (A3) 
G £ (-,0) = Go + e0, (A.4) 

where S,4>, and 0q are given, smooth, time independent functions, and e G (0, 1). It is clear that the 



following lemma, combined with (AT I proves Lemma 3.4 



Lemma A.l. Let e E Uj°l? x n L^H^ be a mean-free global weak solution of ( |A.2 1— ( A.4 1. There exists a 
positive constant C > such that 

||AG e (-,t)|| L2 <C (A.5) 

for all t > 0, and all e G (0, 1), where C = C(||5||^2, ||0||#2, ||0o||#2). 

In [15] we have proven that weak solutions 6 E LfL^, n L^H^ of the unforced AMG K equations ( |1.1[ )- 
(|1.2|>, evolving from merely I? initial data, are in fact Holder continuous for positive time. Moreover, it is 



not hard to modify our argument (De Giorgi iteration scheme) in order to prove the same result if a source 
term S E L p t loc is added, where p > (d + 2) /2. Regarding the higher regularity of solutions, in |[T6l we 
have proven that weak solutions which are C a , with any positive Holder exponent a, are in fact C°° smooth 
for positive time. Clearly if a C°° source term is added to the equations the same result holds. Unfortunately, 
the bound given in JTol for the H 2 norm of the solution may grow with time. The arguments below show 
that in fact uniform in time upper bounds may be obtained. The only e-dependence of the below bounds is 
via ||© e (-, 0)||#2. Since ||0 e (-, 0)||#2 < ||©o||//2 + ||0||/f2, is e-independent, we shall henceforth ignore the 
e dependence of all functions. 
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Proof of Lemma A.l Since ©(-, 0) G C°°, by combining the local existence of smooth solutions with the 
instant regularization proven in lfT5l[l6l . and the uniqueness of smooth solutions, we a priori obtain 

sup ||AG(-,i)|| L 2 < C T (A.6) 

0<t<T 

for all positive time T > 0, where Ct < oo 44> T < oo. We show that Ct may be taken independent of T. 
From the standard energy estimate and the Poincare inequality (0 has zero mean), we obtain 

+ %\m\b < + "livens < \\s\\ l 4&\\li, (a.?) 

which shows that ||0(-, t)\\ L 2 is uniformly bounded in time. Moreover, since S € L^° x , the first part of the 
De Giorgi iteration scheme (cf. El[T5l) proves that 

l|e(-,t)|U» < c||e(-,o)|| L2 (i + t- 3 / 4 ) + c||5|| L oo, (a.8) 



for some universal constant C > 0. The above estimate for t > 1, combined with T = 1 in ( |A.6[ ), shows 
that ||0(-, t) 1 1 l°o is uniformly bounded in time. Estimate ( |A.8[ ) may then be used in the second part of the 
De Giorgi iteration scheme to show that the C a -norm of ©(-, t) is also uniformly bounded in time. The 
latter holds since the Holder norm of ©(-, t) in a parabolic cylinder B r x [to — r 2 , to], only depends on the 
I? and L°° norms of and S on the larger (but finite) parabolic cylinder i?2r x [*o — 3r 2 , to + r 2 ]. 

We have so far shown that ||0(-,i)||cra is uniformly bounded in time, for some a > 0. The boundedness 
of higher order norms is now obtained via energy estimates. Multiplying ( |A.2[ ) by A0 and integrating by 
parts, we obtain 

^H v0 lli2 + < ||«5|| L2 ||A0|| i2 + ||Vt/|| i2 ||V0||| 4 . (A.9) 

To bound the Ty^-norm of 0, we interpolate between C a and H 2 . Using the Besov space characterization 
of the Holder spaces C a = 00 , we claim that the interpolation inequality 

||V0|| L 4 < Clieil^HAeil^ 8 + C||0||^ 8 ||A0|| L / 2 ~ Q/8 ||0||J/ Q 2 (A.IO) 
holds, for any a > and sufficiently smooth 0. Indeed, since we work in three dimensions we have 

||v©|| L 4 < ||5 ve|| L 4 + V ||AjVe|| L 4 



j>0 



< C\\e\\ L 4 + C^2 J '||A i 0|| L / 2 2 ||A i 



l/2||A.QllV2 
L°° 



< cueii^iiAeii^ 8 



+ CY, \\&3®\\% 8 {1 2] \\^M\li) (2 aj '||A i 0|| L o O ) i/2 2-^/ 4 (A.ll) 

where Aj are standard Littlewood-Paley projection operators (cf. [32l and references therein). Upon 
applying the discrete Holder inequality, and using that a > 0, the bound ( |A.11| ) implies estimate ( |A. 10p . 
While the above proof of ( |A.10| ) applies in the case of the domain being M 3 , in order to prove ( |A.10[ ) for T 3 , 
one may use the equivalent Littlewood-Paley description of periodic Besov spaces cf. ||32l Chapter 9] and 
an estimate similar to ( |A.l 1[ ). Since this requires no new ideas, but is notationally heavy, we omit further 
details for the periodic case. 

From ( |A.9[ ) and ( |A.10| ) we obtain the energy estimate 

1 d '|V0||| 2 + k||A0||| 2 < ||5|| L 2||A0|| L 2 + c||e||^ 4 ||Ae||g 4 

+ C||0||^ 4 ||A0|| L r /4 ||0||c", (A.12) 
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for any a > 0. It follows from ( A.12 1, the Cauchy-Schwartz inequality, and the Poincare inequality (since 
O has zero mean) that 



1 d 

2 dr 



|ve||| 2 + 



K 

c 



|VG||| 2 < C K \\S\\ 2 L2 + C K 



MIl 2 



lie 112 



2 Hftll 8 /" 
L 2 H u Hc a • 



The uniform boundedness in time of 1 1 G (• , t) || L i and of 1 1 ( ■ , t) \ \ c a combined with the Gronwall inequality 
proves that ||V0(-, is uniformly bounded in time. A similar energy argument may be carried out by 
multiplying the evolution (jTTTJ with A 2 0, in order to obtain that [| A0(-, £) [| is uniformly bounded in 
time, concluding the proof of the lemma. □ 
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